HIP-2010-12/TH, ITP-UU-10-12, SPIN-10-10 



O 

(N 



X 



Path Integral for Inflationary Perturbations 

Tomislav Prokopec" and Gerasimos Rigopoulos'' 

"'Institute for Theoretical Physics and Spinoza Institute, 
Utrecht University, Leuvenlaan 4, 3584 CE Utrecht, The Netherlands 
^Helsinki Institute of Physics, P.O.Box 64, FIN-00014, University of Helsinki, Finland 

The quantum theory of cosmological perturbations in single field inflation is formulated in terms of 
a path integral. Starting from a canonical formulation, we show how the free propagators can be 
obtained from the well known gauge-invariant quadratic action for scalar and tensor perturbations, 
and determine the interactions to arbitrary order. This approach does not require the explicit 
solution of the energy and momentum constraints, a novel feature which simplifles the determination 
of the interaction vertices. The constraints and the necessary imposition of gauge conditions is 
reflected in the appearance of various commuting and anti-commuting auxiliary fields in the action. 
These auxiliary fields are not propagating physical degrees of freedom but need to be included in 
O ' internal lines and loops in a diagrammatic expansion. To illustrate the formalism we discuss the 

tree-level 3-point and 4-point functions of the infiaton perturbations, reproducing the results already 
' obtained by the methods used in the current literature. Loop calculations are left for future work. 

I. INTRODUCTION 

Perhaps the most remarkable aspect of inflation [l| is its ability to imprint fluctuations on cosmic scales through a 
I confluence of quantum mechanics and general relativity. This connection, first realized more than 25 years ago !2|47| ^ , 
^jj^ has given inflationary theory the impetus which positioned it as the leading paradigm for approaching the physics of 
I— I, the early universe. Since then, the primordial fluctuations have been measured in the CMB Q with ever increasing 
accuracy and resolution 0, [l^ and will be scrutinized even further in the near future [ll| . It is not surprising then 
' that over these past decades a lot of effort has been devoted to fleshing out the predictions inflation makes for these 
fluctuations in a variety of theoretical settings. Since these fluctuations are initially small, of order 10~^ at z ~ 1090, 
QQ the linearized theory of perturbations has been developed to a significant degree and has been used, rather successfully, 
QQ to compare theory to observation. Over the past few years, efforts have intensified to explore inflationary perturbations 
' beyond linear order, mostly in the context of the related non-Gaussianity which has become a significant subfield of 
. cosmological research [l^ . The amount of work that has been done on the subject is by now rather voluminous with 
' many authors examining various aspects. A definitive calculation was performed by Maldacena [13| showing that 
' single field inflation leads to small primordial non-Gaussianities (see also jl^]) with more complicated single- and 
' multi-field models providing more possibilities for larger non-Gaussianity, see e.g. 115.- 311. On the observational side, 
a major effort is under way to develop observational measures of non-gaussianity |10l . |32h41| . 

The foray into non-linear corrections and the associated non-Gaussianity has been motivated by a number of 
reasons. If one neglects the running of the spectral indices, at linear order infiationary theories predict four numbers: 
the amplitude and spectral index of scalar and tensor perturbations, and a variety of different models can coincide on 
these predictions. However, non-Gaussianity can be rather discriminatory for different models due to its much richer, 
and more complicated, structure. As an example, the detection of a significant three-point function would immediately 
rule out single-field inflation, as well as some simple multi-field generalizations, may favor alternative models, or 
provide evidence for the processes involved in heating up the universe after Inflation. Entwined with considerations of 
testing inflationary theory against observations and non-Gaussianity are considerations of theoretical understanding: 
calculating and controlling higher order quantum loop corrections to inflationary predictions, backreaction issues [42| 
and various divergences which appear at higher orders of perturbation theory - see for example 43-56]. 

So far, all the work on inflationary non-linear corrections has been done using the operator language in the interaction 
picture. However, in many branches of modern physics it is the path integral formulation for quantum mechanical 
systems which has proved quite useful. For example, it has been of paramount importance in understanding gauge 
theories and their experimental consequences for particle physics, the theoretical description of condensed matter 
systems and is the preferred language in which modern theories of fundamental physics are formulated and quantized. 



^ The first rigorous and quantitatively accurate treatments of inflationary perturbations were and 0. 



2 



In this spirit, and hoping to throw more hght on the understanding of inflation, here we develop a path integral 
formulation for inflationary perturbations to arbitrary order in interaction terms. An early path-integral formulation 
of linear perturbation theory can be found in [S^l ■ A novel feature of our approach is that there is no need to explicitly 
solve the energy and momentum constraints to a particular perturbative order as is usually done when working in a 
particular gauge. This allows us to obtain the interaction terms, and hence the vertices, to arbitrary order in a closed 
form. 

We start from the more fundamental canonical path integral and show how to obtain the configuration phase space 
path integral, making in the process a connection with the well known gauge invariant linear perturbation theory. 
With a single inflaton there is only one scalar - expressed in terms of the the Sasaki-Mukhanov variable - and two 
tensor degrees of freedom which propagate, as expected, while the vectors completely drop out at quadratic order. 
However, various (real and commuting) auxiliary fields which do not appear in the in-state, the external lines in the 
"in-in diagrammatic" expansion, must be included in computations of N-point functions with iV > 4 or calculations 
involving loops. These fields arise because the constraints are not explicitly solved. Furthermore, anti-commuting 
ghosts arise from the path integral measure which must also be taken into account at a certain loop order. Here we 
focus on standard potential-energy-dominated single-field inflation, but the generalization to more fields and more 
complex theories is in principle straightforward as long as a canonical formulation is available. 

The paper is organized as follows: In section we derive the action for perturbations in a form that includes 
interactions to all orders in a closed expression, discuss the role of constraints and gauge conditions and formulate 
the transition amplitude between two quantum states separated by (background) time t in terms of a path integral. 
In the process the action is written in a form that makes contact with known results from gauge-invariant linear 
perturbation theory [Hs^ and simplifies the computation of the propagators. In section IIIII we discuss the "in-in" 
generating functional which provides the appropriate diagrammatic rules for the computation of the N-point functions 
relevant to cosmology. In section IIVI we apply the above formalism to give expressions for the tree level 3-point and 
4-point functions in terms of propagators, reproducing the results already obtained using an operator formalism and 
a different methodology. We close in section IVl with a discussion of our results. 



II. QUANTIZATION OF INFLATIONARY PERTURBATIONS 



We will consider a single scalar field x minimally coupled to Einstein Gravity. The action in the canonical form [59| 

is (k2 = IGTrGAT = 1) 

J d^x {f'dtg^, + p^dtx ~NH~ NM'} , (1) 
with the hamiltonian H and momentum density 'H* 



W^P^d\~2^p'\ (3) 

where p^^ and are the canonical momenta of the spatial metric gij and scalar field x- They are tensor densities of 
weight one, i.e. p^-' I ^ and Pxl \f9 ^re true tensors. In particular, the covariant derivative in ([3]) is to be understood 
as, VjP^-' = djp^^ +T'^jiP^'' , where F^j is the Levi-Civita connection. Finally, N and iV' are the lapse and shift functions, 
related to the temporal components of the metric tensor as goi = gijN^ and goo = — iV^ + N'^gijN^ , R denotes the 
spatial Ricci scalar, g = det[(7y] and p = gtjP^-' ■ We will be interested in quantizing perturbations around a classical 
inflationary background, so let us set 

P'' - ' (4) 

Px = V^{l+7T^{t,^)), (5) 

9ij = a(t)^((5,y + %(i,x)) , (6) 
X = 0(t) + <^(t,x), (7) 
N = N(t)+n(t,y.). (8) 
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Note that the background vakie of the shift Ni is zero. The action S^^^ for the background dynamics is obtained by 
setting the perturbations to zero 

= J d'xdt {V^ + V^^ NHoit)} , (9) 

where 

The background equations of motion are obtained from S'(°' by varying with respect to a, (p, V and V^, resulting in 

V 

« = -TTT-' (11) 



12a 



(12) 



^--24^ + i^-^"^^' 

= -a^V^, (14) 

where V — V{4){t)), V^) — dV{(j))/dip and we introduced the 'dotted' derivative, d = N^^da/dt, etc. Furthermore, 
variation with respect to N gives the constraint Holt) = 0, or equivalently. 



24a 2a' 



(15) 



The well known equations of flat Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology are obtained with the 
identification 

Indeed, inserting the solutions for the momenta V and 7^^ into (IT5] - [T5|) one obtains, 

(17) 

(18) 
(19) 

(P0 = Y - ^) ■ (20) 

Recalling that a dot refers to d/ (Ndt), we see that these equations are time reparametrization invariant. For example, 
their form in cosmological and conformal time is obtained simply by choosing the lapse iV = 1 and N — a, respectively. 
Of course, only two out of these four equations are independent. Indeed, from the first equation and any one of the 
other three, one can derive the remaining two equations. To recover the dependence on the Newton constant Gjv one 
should reinsert — WttGn = 1 into the r.h.s. of the first two background equations (fT7|) and ((T8|). 

In such a model, inflation takes place when the fleld is rolling slowly on the slope of its potential. We therefore 
recall the 'slow roll' parameters, which for the purpose of this paper we deflne as, 

*^="ff2 = o— ' ^=^17l- (21) 
2 p0 H(j) 

Notice that in the slow roll approximation when e ^ 1 and 77 1, e and rj reduce to the standard slow roll definitions, 
e — > (Vcfi/V)^ and 77 {Vift'/^) ~ ^- For the rest of the paper we wfll be using H instead of V and V,p instead of (j). 
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A. The action for perturbations 

We can now proceed to obtain the action for the perturbations by making the replacements (|4]-[8l) and expanding ([TJ . 
We ignore terms hnear in perturbations since they are multiphed by the background equations of motion (jl7H20l) . 
which we assume to hold. We can therefore write the complete action for the perturbations in canonical form as 

5port = J d^xdt (v^TT^dt^ - 2a^Hn'^dtKj - Hpcrt + nCo + N,C,^ = Sp + Si (22) 

5p . /Ad,{p.,,a,„-2„W^4,..,-«p + nCi'> + iV.C«} (23) 

Si = j d^xdt { - -Hi + nC}^ + N,Cf^^ , (24) 

which completely specifies the dynamics of perturbations to all orders in the interaction terms. For later convenience 
we have split the full action iSpert ((22|) into the free (quadratic) part Sf P3I) and the interactions Si ((24)) . which include 

cubic, quartic and other higher order terms in perturbations. Cq — c'^^ + C^^ = — "H^^^ /N — Hpcrt/N stands for the 
hamiltonian constraint, which is split into the linear part (denoted by the superscript (1)) and quadratic and higher 
order parts (denoted by the superscript > 2 and the subscript pert). We have also split the momentum constraint as, 
Ci = Cj^^'' + Cf^ = — "H*. Analogously, "Hpcrt = 'Hp + "Hi is the full hamiltonian dictating the time evolution of the 
perturbations, split in Eqs. (j23H24p into a free (quadratic) and an interacting part (qubic, quartic, etc). The Poisson 
bracket algebra of the system of perturbations can be calculated to close 

{Hpcrt(x), Co{y)} = - (a(x)af - a{y)d^) <5(x - y) + ^J(x - y) , (25) 

{Hpert(x), C.(y)} = Co(y)9,r5(x - y) + ^<5(x - y) , (26) 

{Co(x), Co(y)} - (Q(x)af - aiy)dn ^(x - y) , (27) 
{C.(x), Co(y)} = Co(x)9f<5(x-y), (28) 
{Q(x), C,{y)} = {C.iKjaf - Q(x)9j') 5(x - y) , (29) 

where the fundamental Poisson brackets are seen from (|22p to be 
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My,t), ^^(x,i)} = — 5(x-y), (30) 



{hki{y,t),7r'^i^,t)} = --l-{S,kSji+SuSjk)5i^-y). (31) 

We see that from a canonical point of view, inflationary perturbations form a constrained Hamiltonian system [6C 

What is the number of dynamical degrees of freedom in such a system? Following the point of view in [61|, we 
see that the seven fields {ip, hij) appearing in (|22p are not all independent degrees of freedom: n and Ni are lagrange 
multipliers without dynamics of their own, imposing under variation the conditions Cq = and Cj = which 
constrain the evolution of the system to take place in a lower dimensional hypersurface of the full 14-dimensional 
phase space. The functions n and Ni are not determined by the dynamics and are arbitrary; in this particular case 
they parameterize the freedom in choosing how to break spacetime into spatial hypersurfaces and time. Thus, the 
dynamical equations contain four completely free functions which need to be fixed by the imposition of four further 
gauge conditions: Qa{hij,ip,'K^-' — 0. These gauge conditions can be chosen at will, up to the requirement 

{Qa, Q0}=O, Det{Q„,C/j}^0. (32) 

We see that the physical phase space has dimension 14—4—4 = 6. Indeed, as is well known and will be re-derived below, 
there is only one dynamical scalar and a transverse traceless tensor propagating in single field infiation, corresponding 
to 3 degrees of freedom and a 6 dimensional dynamical phase space. Note that once a choice of Qa has been made, 
the lagrange multipliers Na = {n, Ni} can be determined by the imposition of the consistency relation 

Qa = {Hpcrt, Qa} " Np{Cp, Qa} « , (33) 



For an early treatment of inflationary perturbations as a constrained Hamiltonian system see [62 
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where the symbol w imphes that the constraints = are imposed after the poisson brackets are evaluated. 
Condition ([5^ then insures that equation ([55)) can be solved for N^. 

Since the number of dynamical fields in the system is much less than the apparent dimension of its phase space, 
it proves convenient to explicitly separate out the 3 propagating degrees of freedom from the rest. We do this in 
the next paragraph, leading to equation (|45|) . Let us note that the physical degrees of freedom could be isolated by 
actually solving the constraints, thus expressing some of the canonical variables in terms of the others, and pl ugg ing 
the solution back in the action. This is feasible for the free theory, where the interactions Si in are ignored pT\ . If 
the interactions are included however, the solution of the constraints can only be found by an iteration procedure and 
leads to rather cumbersome expressions. The approach we follow in this paper dispenses with the need to explicitly 
solve the constraints, although by direct comparison it is equivalent to such a solution at the linear level. 

Recall now Eqs. (|22HM)) . where the action for perturbations was split into the free and interacting parts. Explicating 
the free action Sp the free hamiltonian reads. 



•p2 



.3 r 



2a3 



hn, 



N- 



hllhl + Na( -hV^h ~ -hd,d,K, 
4 V 4 2 ■' 



-hijdidihji 



(34) 



where we have used the notation Aijki — 5ik5ji + 5ii5ik — SaSki, and 



'JikOjl 



'JilOjk 



'JijOkl, 



(35) 



C 



(1) 



'-d^^~AaH(dnT'' + dihi, 



-d,h 



(36) 



and where = Sijdidj. From now on a superscript/subscript > n will indicate that only terms of order > n 
in the perturbations are included, such that Cj^^ denotes contributions that are quadratic or higher order in the 
perturbations. 

Next, it will be convenient to use the standard scalar- vector-tensor decomposition of the spatial metric perturbation. 



y 3 



, TT 



with 



1 



0. 



Furthermore, we will decompose the shift into the longitudinal and transverse components as, 

Ni^diS + NT, with diNT^Q. 



(37) 

(38) 
(39) 



After a series of manipulations, which we present in appendix B, the free action S-p can be diagonalized. In summary, 
if we define 



(40) 



and perform the following linear shifts in all the perturbation fields: 

n = n + Lq 



u 

-U 

L. 



(41) 
(42) 
(43) 
(44) 



where Lg, Li, and are given in Eqs. ()139H143| of appendix B, we find that, up to boundary terms, the free 
action takes the form 



y (fxNdta^< 




r.2 fdiw\-^i 
" [—) . 


1 




.18 




+ 4 



-2Hhpl-2H^p'^A,jkip'" 



2(3 - e)H' 

iV2 



Njy^Nj -^^{^^Sf^Y (45) 
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In this form, the action clearly shows that the perturbations which propagate are one scalar degree of freedom 
w and the 2 degrees of freedom of the transverse traceless tensor hj^^ . All the other fields are auxiliary fields 
without dynamics of their own. The first line of (|45l) is well known from the gauge invariant theory of cosmological 
perturbations [ssj ^. The fields w and hj^^ in Eq. (|45p are gauge invariant under linear gauge transformations and 

time independent on long wavelengths It is interesting to note that all the auxiliary fields, both h, Nf and S as 
well as the shifted momenta p*-' and p^p are also invariant under such linearized transformations - see Eq. (jl47p in 
appendix B and the discussion that precedes it. 

The variation of the action ps)) w.r.t. tt^, tt*^ , n, Nj and S yields the on-shell relations, p^p — 0, p'-' = 0, ri = 0, 

NJ = and S — 0. The relations p^p = Q and p*^ — translate into the relation between the momenta tt*^ and 'n^p 

and the time derivatives of hij and ip. On the other hand, n = 0, V'^Nj — and V^S* — translate into the solutions 
of the linear constraints. From Eqs. () 13911140]) and ()143p we see that the lapse perturbation n and the shift Ni are 
then given in terms of the fields (iy9, hij) and their one time derivative or one or two spatial derivatives. Thus, the free 
action (1451) contains all the elements of the well known gauge invariant treatment of linear inflationary perturbations 
and the transition between the hamiltonian and lagrangian formulation of the problem. 
Interactions are encoded in Si given in Eq. (|24|) . with 

Co-' = -i(HF+Hi) (46) 

and 

= - '^{f)>idiV - '^^^fdi^ - AaH \^{f)>i(d,h,i - ]^dih) + f(d,hir - l^^/i^,) ^^''j , (47) 

where T-L-p and "Hi are given in Eq. and Eq. (I159|) of appendix D, respectively. After all the variables are expressed 
in terms of the redefined fields, n, Ni, p^^ and p^p it becomes clear that Si cannot be written solely in terms of w, hj^ 
and the auxiliary fields which we will collectively call a. Indeed, if is exchanged for w in Si, ip ^ ^{h — V^/i — w), 
one obtains an interaction action ([M|) which is a functional not only of w, hj^ and cr, but also of h, h and hj (for 
a more explicit form of S\ see Eqs. (|159|) and (^Sl - HTl) '). Thus, the interaction terms seem to contain four more fields 
than those appearing in Sp (|45p which determines the free dynamics. This is of course related to the aforementioned 
arbitrariness in the time evolution of the system. The imposition of the four gauge conditions Qa = cures this 
arbitrariness by expressing h, h and Kf in terms of the fields appearing in Sy- This observation also demonstrates 
the necessity to consider non-linearized gauge transformations when interactions are included. Notice that Sp + Si 
is explicitly invariant under linearized transformations, and hence it would seem that the action as a whole is not 
since Si explicitly depends on the gauge chosen, even at the linear level. This contradicts the well known exact 
gauge invariance of the action^. The resolution of this seeming paradox lies in the observation that the exact gauge 
invariance of the action implies that for interaction terms of order up to n, one must he able to absorb any such 
apparently gauge dependent terms by non- linear field redefinitions: gauge transformations of order n — 1. An explicit 
example for cubic interactions has been worked out in [13]. With this observation, quantities that have been calculated 
in different gauges can be compared using the appropriate field redefinitions. 

B. The path integral 

The results of the previous section can be used to obtain a convenient quantization scheme for the perturbations, 
including the interactions. According to the above discussion, inflationary perturbations are described by a constrained 
hamiltonian system with action (j22p . In general, the quantization of such systems is easily formulated in terms of a 
path integral [6l|. Indeed, the transition amplitude from state |-4;iy!i) sA t — t^ to \B]tB) at time ts can be written 



^ The scalar perturbation is usually expressed in terms of the Sasaki-Mukhanov variable v = a(Lp — z(h — V^/i)) with z = T'^/lda'^ H], and 
the corresponding free action, when written in conformal time, is shown in Eq. 11148 I I of appendix B. 

If the gauge ip = is chosen, w represents the curvature perturbation and is also constant during inflation on long wavelengths to all 
orders in perturbation theory 
^ Up to boundary terms which we ignore in this paper. 
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(/l>j,¥')(tB)=B 

{B;tB\A;tA) ^ J [Vh,,VipV7r''V7r^VnVN,]Y[S[Qc.]\Det{Q^,Cp}\e'^-^'-, (48) 

{hij,ip){tA)=A 

where 5port is given in (|22l - [24l) . integrated over time t' e [t^jtg], and we have schematicahy denoted the initial (A) 
and final (B) field configurations by the limits of the path integral. The transition amplitude (pS)) is the fundamental 
quantity for calculating quantum correlators of any kind and the starting point for calculating the correct lagrangian 
path integral. According to [61], it is independent of the choice of gauge conditions Qa, as long as the boundary states 
are appropriately defined in an invariant manner. We shall address this issue in a forthcoming publication. Note the 
appearance of a non-trivial measure containing |Det {Qa, C;?}!, the determinant of the Poisson brackets between the 
constraints Ca and the gauge conditions Qa- Equation p8)) is a self-contained starting point for doing perturbative 
calculations. By splitting the action as above in 



5pert = 5f + 5i , (49) 

it is possible to evaluate the path integral (|48p perturbatively, expanding in powers of the interaction action, after 
including an appropriate representation of the determinant. The free action Sf then determines the propagators while 
Si gives the various interaction vertices, completely defining a diagrammatic expansion. 

The focal point of the previous section has been the free action and (IMH55|) that can be brought to the diagonal 
form (l45t from which propagators can be easily obtained. The field redefinitions (|4TH44)) are linear shifts in all the 
fields of the action which do not affect the form of the functional measure in the path integral. More precisely, they 
contribute time dependent - but field independent - factors which cancel in our later computations. The transition 
amplitude can therefore also be written in terms of the redefined fields as 

{B;tB\A;tA) ^ J [VwVhf^VaVhVhVhJ]Y[S[Qa,] \Det{Qa.,Cp}\ e'^^^'"'''^^^'''^+'^\ (50) 
(/if/,«,)(u)=^ 

making it evident that the physical propagating degrees of freedom of the system are the quanta of w and hfj^. 

Va = Vp^^Vp^VfiDNf'DS collectively denotes integration over the 11 auxiliary fields. One can note that the full 
action is quadratic in pij and which could in principle be integrated out. However, it proves convenient to keep 
them along with n, Nf and S. The "propagation" of all these auxiliary fields is determined by the second line of p5)) 
and is trivial. We here use the term propagation for the auxiliary fields rather loosely since they are not attached to 
wave operators. Since w and hj^^ are the dynamical degrees of freedom, it is clear that the initial and final states 
should be determined in terms of their free quanta. 

To complete the derivation of the more workable expression for the path integral, we are left with the task of 
evaluating the determinant |Det {Q^, C;3}| in ([50)) . Denoting 

r!a/J = {Qa,C0}|g^^O, (51) 

where the gauge condition — has been imposed after evaluating the Poisson brackets, we can represent the 
determinant as a path integral over anticommuting scalar fields - Faddeev-Popov ghosts - t^q, and rja ^ 

\Detn^p\ = j [Vfj^Vrj^] e'/ . (52) 

The transition amplitude is then written as 

{B;tB\A;tA) ^ J [VwVhJ^VaVhJVhVhVfjVrj] ]j6[Qc,] e'^^ + '^'^''-^^-^'^^ + '^' + '^ ^ (53) 

{h'[/;w){tB)=A 

where we have split rtap = ^Q^ + ^^i^ into a free part, determining the "propagation" of ?7q, and r/a, and an interaction 
part which contains the other fields and couples them to fja and r]a . Two common gauge choices Qa and the resulting 
matrices fla/B are discussed in the appendix. We now have a complete description of the transition amplitude amenable 
to perturbative calculations. Equation (j53p is the central result of this paper. 



See for example Eqs. (5.1.31-32) in Ref. [63 
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III. THE FUNCTIONAL IN-IN FORMALISM AND EXPECTATION VALUES 



The previous section described the derivation of a path integral representation for the transition amphtude between 
states of inflationary perturbations in the past and the future, separated by time t. The perturbative evaluation of 
such amplitudes between the far past and the far future leads to the usual Feynman diagrams for the calculation of 
S-matrix elements, i.e. scattering amplitudes between in-states and out-states, useful in deriving predictions for the 
scattering experiments of particle physics. However, in cosmological applications the quantities of interest are not 
scattering amplitudes between predefined states (a boundary value problem), but expectation values of operators at 
a time t, given an initial state |in,t^) at time (an initial value problem). For the purpose of this work, this state 
is defined in terms of the freely propagating modes of w and hj^^ at very early times: — ^ i-oo, and is taken to be 
an appropriately defined vacuum state |r2,<_oo). More generally, one has to begin the evolution from a finite time 
and with an initial state chosen such that it includes corrections due to the interactions, such that all field correlators 
are finite at the initial hypersurface. In this setting, the calculation of expectation values leads to a modification of 
the standard diagrammatic rules of QFT. Here, we sketch the functional in - in formalism for obtaining expectation 
values and arrive at the appropriate diagrammatic rules. For more details see [g^, HI] as well as the appendix of [4j]. 



A. The in-in generating functional 

In analogy with standard QFT, one defines the in-in generating functional Z[J_, J+] which contains two arbitrary 
and independent sources J-|_ and J_ 

Z[J_, J+] = ^^{^,t-oo\a,tout)j.{a,tout\^,t-oo)j+ , (54) 

a 

where {a,tout\^,t-oa)j denotes the transition amplitude from the aforementioned vacuum state at time <_oo to a 
state |a,io«t) at time tout in the presence of a source J. Expression ([5l)) can be interpreted as going forward in time 
in the presence of a source J+ and then returning back in time in the presence of a source J_ . This explain the term 
"closed-time-path formulation" that is sometimes also used. However, in what follows such an interpretation is not 
needed. The final time tout can be taken to be any time later than all the times of interest in the calculation; here we 
take tout to be in the far future, tout +00. Obviously, Z[J, J] — 1. 

The transition amplitudes appearing in (j54p can be expressed in terms of path integrals. Collectively denoting the 
relevant fields by ^ we have that 

Z[J_, J+] = J [p^^p^_]e*SK+]-.SK-] + u+5+-»j-C- ^ (55) 

with boundary conditions appropriate for the state jfi, i_oo) at i_oo and the constraint that S,+ {tout) — £,-itout)- Once 
Z[J^, ,/+] has been calculated, expectation values can be obtained by taking variational derivatives: 

Z[J-,J+] (56) 



^iSJ-{yi) ~i6J-{ym) i6J+{xi) i6J+{xn) 



J-l- = J_=0 

= {Q; t_oo |{T^(yi) • • • ^(2/™)}{rC(xi) • • • a^n)m, t^^) , 



where T denotes operator time-ordering, T^(x)^(x') = 9{x'^ — x"^)^{x)^{x') + 9{x"^ — x'^)^{x)£^{x'), and T operator 
anti-time-ordering, f^{x)({x') = e{x'^ -x^)^{x)^{x') + e{x° -x'°)^ix)clx'). 

The exact calculation of Z[J+, J_] is impossible for interacting theories but it can be evaluated perturbatively. 
Writing the action as a free part plus interactions 

S[^]= J d'x^^Vo^ + Sii^], (57) 
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we have ^ 



11 d4:,dVi(a„-u_)(x)(:t:;:t:)(a;;a;')(_^,;S"_)(^')} , (58) 

where Vq denotes the operator acting on ^ as imphed by the free action (j45p . The last equahty is derived by shifting 
the fields 

(t)-(t)-(:A!::^!:)(-fL), (59) 

where the Keldysh 2x2 matrix of propagators is defined as the matrix inverse of Vq, 



(^) A + A ]{x■,x')^^S\x^x'). (60) 



Vo \ ( iA++ iA+ 
-Vq J ^"^f \ zA_+ 2A__ 

The J-independent path integrals in (|58|) contribute phase factors which cancel. The diagonal elements in (|60|) are the 
(Wightman) Green functions of the operator Dp which appear in usual S-matrix calculations in QFT, while the off- 
diagonal elements are solutions of the homogeneous equations, VqiA |_ = VqiA^ = 0. These off-diagonal elements 

in (j60p are necessary in order to impose the boundary condition ^-|_ = ^_ at time tout on the J-dependent parts of ^+ 
and ^_ . 

The explicit form of the A propagators can be found by using ([55)) and ([SS)) . after setting Si — > 0. We see that for 
propagating fields the appropriate elements of the Keldysh (bosonic) propagators in (1601) are 

iA++{x;x') = {^\Ti{x)£,{x')p) ^e{x'^-x'^)iA^+{x;x') + e{x"^-x'^)iA+^{x;x') (61) 
iA^-{x;x') = {Vt\fi{x)(,{x')\n)^9{x°-x'^)iA+-{x-x') + e{x'°-x°)iA^+{x;x'), (62) 

where the Wightman functions iA_| and iA |_ are defined as, 

iA+_{x-x') = m{x')axm, iA^+{x-x') = {n\£,{x)i{x'm, (63) 

and denotes the physical vacuum state If the fields are non-propagating, the Wightman functions vanish, 
A+_ = = A_+ and Eq. ^ is easily inverted, iA±± = ±{Vq)-'^ [x)i5'^{x - x'). 

B. Propagators 

Let us now focus on the problem at hand and in particular on the propagators that correspond to Sy given in (j45p . 
The propagators for w 

iA1+{x\x') = {n\Tw{x)w{x')\9) , iA'^_{x-x') = {n\fw{x)w{x')\^l) 

iA'^^{x]x') = {n\w{x)w{x')\n) , iA1_{x\x')^ {Vl\w{x')w{x)\n) , (64) 

and for hf^ 

zAf^ix;x') = {n\Thj;;'{x)hin^'m), ^^l'!'{x■,x') ^ {n\fhj;'{x)hlf{x')\n) 

.A!^"t'(x;a;') = {n\hj;;'ix)hl^{x')\n), zAf J {x; x') = {n\hj;;' {x')hl^ {x)\n) (65) 
can be obtained by solving the equations, 

(-StZ^flS^ _,_ iVz^aV^) lA^^ix; x') = ±iS^ix - x') , (^-dtz^a^^^ + Nz^aV^^ jA^^(a;; x')^0, (66) 



Formula II58I I applies also for the (fcrmionic) ghost fields, as can be seen e.g. from Eq. (5.1.32) in Ref. [63l| . where iA sue the ghost 
propagators. 

* In statistical field theory pure states are replaced by more general mixed states (see e.g. [6Q . 67]), and the evolution of a system is generated 
by the density operator p{t). In this case the Green functions are defined by the more general relation, iA^^ {x;x') = Tr[p(t)T^(x)^(x')], 
etc, where the definitions must be picture independent. For example, if the density operator p{t) is in the Schrodinger picture, the fields 
^{x) must be in the Heisenberg picture. 
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and ^ 

i-dta^^+NaV') lAf^ix; x') = ±{P,kPji+PuP,k~P^jPki)i5\x-x') , {-d,a^^+Na^^) i^l'^{x; x') = , (67) 

where z = v^/3, Pij = Sij — didjfV^ denotes the transverse projector, and the operator is defined by Vx^(5(x — 
y) = ■iAo(x; y) = — l/[47r||x— y||]. The projectors on the right hand side of (|571) project out the potentially unphysical 
(scalar and vector) 'degrees of freedom,' which are not transverse or which are not traceless, assuring the correct 
counting of the physical degrees of freedom. Due to the infrared problems which plague Eqs. (|66H67I) in accelerating 
space-times, the corresponding solutions are not simple and need to be handled with due care [52, 53, 68J. Here we 
just note that, when the fields in (j64[) and (|65p are Fourier transformed to the spatial momentum space, 

g-.k.x ^ ef^{k)[hit,k)a^^^ + h*{t,k)a+_^] (68) 

a=+,x 

and written in conformal time (N o,{ri),t — >■ r/), Eqs. (I66H67|) simplify to, 

{-Oydri - a^k^) i^i'iim V'; k) = ±{P^kPol + PaPjk - PrjPki)iSiv - v') , 4 = % - ^ • (69) 

We arrived at these results by choosing the physical state that is destroyed by the annihilation operators Ok and 
<iQ,k in (1551) . flkiri) = 0, tta^kl^) — 0. The solutions of Eqs. can be expressed in terms of the mode functions 
10(77, k) and h{r], k) as, 

i^l^ ^ 9{ri - r]')iA'^^ + 0{r]' - r])iAl_ , iA'^_ ^ w*{ri,k)w{r)' ,k) , iA!^^ = w(r/, fc, )u;* (r/, fc) 
zA!f+" = % - ryOiA'^r + ^(^' - '?)«^+-"" 

^^n.n ^ J2 e'^{k)eUk)h*{i^, k)h{i, k) , ^A'ir = ^ e%{k)eU'k)h{r^. k)h*{rj', k) , 

a a. 

where w{r]^ k) and h{rj^ k) satisfy, 

(-52-k2 + ^^)(azu;(7y,fc)) ^0, (-92_k2 + ^)(a/i(ry,fc)) =0, (70) 

where efj denotes the graviton polarization tensor. ^'^ These equations can be solved for certain choices of the scale 
factor a — a{ri). For example, when e = —H/H^ — constant, a = [—{l — e)HQrjY/^^~'^\ where Hq — H{rio), such that 
a" /a = (2 — e)/[(l — e)^?7^] and (|7D|) can be solved in terms of the Hankel functions, 

w{k,v) = -J^\a^{k)Hl^\k\r^\)+/3Uk)Hi'\k\v\)\, ^=7^^,^ |a^(fc)P - |/3^(fc)|2 = 1 
az V 4 L J 2(1 — e) 

KKv) - -A/^prU(fc)ilW(fc|7y|)+/3,(fc)i7(2)(fc|7y|)l , K(fc)p - |/3,(fc)|2 ^ 1 , (71) 
a V 2 L J 

and where ^-u, — Vh = v = (3-e)/[2(l— e)], and the coefficients {aw{k), I3w{k)) and {ah{k), /3h{k)) of positive and negative 
frequency mode functions are chosen such that the real space propagators are infrared finite pa [s^ . Alternatively, 



The normalization on the right hand side can be checked by noting that, if the graviton were not transverse or traceless, i.e. just an 
ordinary 3x3 symmetric matrix, then the propagator equation would read 

Eq. II67I I is obtained from this equation by removing the longitudinal and the trace degrees of freedom. 

The polarization tensor ef- satisfies, Ec=+,x ^g-^H = (l/2)[AfcPj! + PuPjk - P^JPkl] and ^fj^j = ^"'^ ■ 
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one can place the Universe in a large comoving box [43|, |45|, |52|, |68| , or in universe with spatially positively curved 
sections, which effectively remove the far infrared modes, thus rendering the state infrared finite. For non-constant e, 
solutions can be obtained in the slow-roll approximation, by using the de Sitter mode functions around the time of 
first Hubble crossing [l^, or the methods of |69| . 

The two-point correlators for all of the auxiliary fields, p'-', p^, fi^Nf and S, can be straightforwardly obtained. 
It is clear that the off-diagonal elements (Wightman functions) of (|60l) vanish, while the diagonal elements are either 
proportional to delta functions or involve spatially non-local operators. In particular, the Green functions for the 
auxiliary fields as implied by (j45l) are 



1 



1 

N 



tS\x~x') 

(SikSji + 5ii5jk - 25ij5ki)iS'*'{x - x) 
i5Hx-x') 



4(3 - e)a3i/2 
TaNV-^PijiS'^ix-x') 

±laN(i-e)V-^i5\x-x'), 



(72) 
(73) 

(74) 
(75) 
(76) 



with all the zA^ = for the auxiliary fields. 

Let us now turn to the ghost fields. Their propagators satisfy the equations of motion of the general form, 

^Ip (*A^±±) {x- x') = ±z5^,5\x - x') , [lAl^]. {x; x') = . 



(77) 



Whether the ghosts are propagating or not depends of course on the choice of gauge conditions and in particular 
whether the latter contain time derivatives. In the former case we have 

{iAl^)++{x;x') = (f^|r77„(x)r;^(x')|f^) =0(a;°-x'")(zA^^)_+(a:;a:') + 0(x'"-a:°)(zA'^^)+_(:r;a;') 
{^Al^)..{x■,x') = (f^|f77„(x)^^(x')|f^)=0(a;°-:r'°)(zA^^)+_(a:;a:') + e(:r'°-:r:°)(zA';i^)_+(:r;a:'), (78) 



where 



{zAl^)+.ix;x') = - {n\fjp{x')r^^{x)\n) , (zA^^)_+(x; x') = {n\7^^{x)np{x'm 



(79) 



Note that because they are anticommuting, their corresponding time ordered and anti-time ordered propagators differ 
slightly from the bosonic counterparts ([6n - |62|) . Since common gauge choices do not lead to propagating ghosts, the 
corresponding propagators have simpler form than the one given in (j78H79p . Consider, for example, the tensor gauge 



(80) 



Imposing the gauge conditions Qa — sets h ^ h ^ hj — and hence the metric perturbation contains only the 
(transverse traceless) tensor part of the spatial metric, while w represents the perturbation of the infiaton. In linear 
theory it corresponds to the uniform spatial expansion gauge, or the comoving gauge. 

From Eq. (|156p in Appendix C it follows that the free part of the ghost operator fl^^ is non-dynamical in the tensor 
gauge; it reads: 



-6H -2a-^d] \ 

-a-2(5,^.v^ + iafa;) 



(81) 



Eq. l|73|l is obtained by noting that the inverse of A^ji^i is (A~^)iji,i = (l/4)((5j(;5j; +<5ii<5_,7; — 2&ij6i^i) in the sense that Ayi,; (A~^)f^i^„ = 
i((5im(5j„ -I- SinSjm). That this indeed is the correct definition of the inverse can be checked by calculating, for example, the propagator 
for the (scalar) trace p = Sijp''^ of p*-' , for which the free action I I45II yields J d'^x[2Na^ /3]. The corresponding propagator is, of 
course, iA'^_j_{x; x') = ±[3/ {4:N H'^)]iS'^ {x — x'), which agrees with the double trace of Eq. II73I I. obtained by applying SijSj^i. One can 
analogously check that the propagators for the other components of p*^ are correct. 
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The ghost equations, when written in a condensed matrix notation: • {iA^)±± = ±i5^, f]^ • (iA'')±zp = 0, imply, 



and (iA'')±zp(x;x') = 0. 

The second gauge we consider here is the uniform field gauge, 

Qo = ^, Q^ = ^J {h^j-^h^ . (83) 

In this gauge = fixes = h = hj ~ and hence the metric perturbation contains the tensor and the spatial trace 
perturbation. This means that w now represents the perturbation of the local spatial volume, or the spatial curvature 
perturbation. From Eq. (jl58l) we see that, also in the uniform field gauge, the free ghost operator is nondynamical, 

-2^H 



implying the following 2-point ghost correlators 

{iA^)±±ix;x')^(j^'^'" ^^^^^ 1 ]t6\x-x'); (iA'')±^(a;; a;') = . (85) 



4 vi y V 



Notice that the gauge (155)) is not suitable when one is interested in the de Sitter limit e — 0, since the ghost 
propagator ([85]) appears singular in that limit. That problem can be fixed by simply replacing Qq in ([SS)) by Qg = ip/4>. 
In order to make the whole action (P^VE^ regular, one also needs to redefine the canonical momentum of (p as 
VipTTip — ^ TTip, and similarly —{^/3)w — ?> ip, such that the corresponding propagators ([72)) and (l66l) become regular. 
When these changes are exacted, the de Sitter limit of our path integral ([53]) will be regular in the uniform field gauge. 



C. Diagrammatic rules 

After expanding (|58p in powers of Si, we are led to a diagrammatic expansion for Z[J+, J_], analogous to that 
of the vacuum-to-vacuum amplitude in the presence of a source J of standard QFT. Diagrams are now composed of 
external lines coupling + and — currents at external points to + and — vertices, as well as internal lines coupling -|- 
and — vertices with each other. Each vertex carries a +i or ~i factor respectively, along with the relevant interaction 

factors. Vertices and currents of the same valence, {+, + and — , — ) are joined by lines corresponding to A++ and A 

propagators, respectively, while those of opposite valence (-I-, — and — , +) are joined by A^i and A ^ propagators, 

respectively. Finally, one integrates over the coordinates of internal points. 

It is important to keep in mind that, as mentioned above, the 2-point correlators of the auxiliary fields and of the 

ghosts in common gauges satisfy A_| = A |_ = while the A_|_+ and A are either delta functions or involve 

non-local spatial operators of the type (5(i — t')/||x — x'|| or even (5(i — t')V^^||x — x'||~-'^ = —{A7r)~^d{t — t') J cfiz\\x — 
z||^^||z — x'||~^ as in Eqs. (|75H7S)) . Since there is no actual propagation in spacetime, these correlators always 
correspond to effective vertices, either local or non-local. 

Once Z[J-,J+] has been expressed as a sum of diagrams to the desired order in the interaction or in loops, 
variational derivatives will then give expectation values according to . Note that in case all fields in ([5B)) are taken 
at the same time, the valence of the currents is immaterial; here we will be using J+, so that e.g. 

mit,^)at,ym) = / / Z[J^,J+] , (86) 

idJ+(t,x) idJ+{t,y) j-=j+=o 

and 



i(SJ+(t,x) iSJ+{t,y) i5J+{t,z) 



All of the above mirrors the corresponding discussion for standard QFT transition amplitudes. We can thus 
summarize the diagrammatic rules for the calculation of an n-point function (51|^(t, x)^(i, y) . . .^{t,z)\fl) as follows 
(see the appendix of for a complete definition of such rules): 
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• draw the usual Feynman Diagrams with all external points having a + valence; 

• all internal vertices are either + or - and carry a +i or —i factor, respectively, along with all the relevant 
interaction operators; 

• connect external points with vertices and vertices with each other by the propagators of the appropriate valence; 

• integrate over the temporal and spatial coordinates of all the vertices. 



D. The tensor gauge and cosmological perturbations 



The N-point functions of the trace of the spatial metric Xi)/i(i, X2) . . . Xjv)|f2) on comoving spacelike 

hypersurfaces are significant cosmological observables, particularly in relation with the non-Gaussianity of cosmological 
perturbations. During inflation, such quantities can be directly computed in the uniform field gauge, since in that 
gauge w ^ h. However, the interaction terms in Si are simplest in the tensor gauge and the question arises of 
how to relate results in these two different gauges. The issue is that, although the transition amplitude ([53| is 
independent of the choice of gauge conditions, simply coupling w to a current and calculating expectation values does 
not automatically lead to a gauge invariant expression. This is obvious since w is not invariant under non-linear 
transformations and adding a Jw term in the action explicitly breaks its gauge invariance. A more physical way 
to understand this is to note that a time slice defined by some fixed value t of the background time corresponds 
to physically different hypersurfaces in different gauges. Thus, one cannot simply compare {w^ (t)) for fixed t in 
different gauges because it refers to physically different quantities evaluated on different timelike hypersurfaces. Such 
difficulties can be circumvented however if we define the observables of interest and the timelike hypersurfaces on 
which they are to be calculated in a geometrical fashion [70j . 

So, let us choose the tensor gauge to obtain the interaction terms. The physical interpretation of w is that it 
represents the quanta of the scalar field: w = —{3/y/e)ip, while the metric perturbations are transverse traceless 
tensors. Interactions between these degrees of freedom are specified by Si in this gauge, given explicitly in section ITVl 
below. The transition amplitude that enters in (|54|) is then written as 



where the operator 17^^ is the one appropriate for the particular gauge choice - see ([5T|) and (|156l) - and the vacuum 
states have been obtained by considering large time intervals and the appropriate ie prescription. The expectation 
value of any product of w's can now be calculated using the diagrammatic rules described above, remembering that 
the meaning of correlators in this gauge is 

3. A. 



(f7|ii;(i,xi)ii;(i,X2)...«;(i,XAr)|fi) = (^-^j x)(^(t, xa) . . . XAr)|f7) . (89) 

Let us now consider the hypersurfaces defined by the condition if = Q, which are of course different from the 
t = constant hypersurfaces in the tensor gauge we are using. We are interested in the determinant Det[7y]ip=o of the 
3-metric ^ij induced on these hypersurfaces which measures their local expansion. Writing 

Det(7y)^=o = a{tfe^^ (90) 
we can calculate the expectation value of any N-point function of Q from 

(»|C(^,xi)C(^,X2)...C(^,x^)|») ^Cf ... O ^f i ^ ^ ■■■^{ ■,1 I ^ ^ )z[J^^J+\ , , '(91) 

\lOJ+(t,Xi)/ \ldJ+[t,X2) J \ldJ+[t,XN) J J_ = J^=Q 

where on the r.h.s. C(''^) h^-s been expressed as the appropriate function of w and all occurrences of w have been 
replaced by 5/[i5J+{t, x)]. For example, to second order in perturbations, C, is related to w in the tensor gauge by (c/. 



e.g. [13|): 



36 C ~ 6w + ye — -77J up' 

1 . 1 1 f,r. s,^ s didi/,^ ,\\ . ^ di .\ e didi ( ,^ ^ / c?, 

-\-—ww — 7- 



-^h~^a.ai^. (92) 
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When restricted to long wavelengths, formula ([9T|) could be considered as a quantum mechanical generalization of the 
AN formula. Analogous expressions can be used for relating any quantity of interest to results obtained in the tensor 
gauge. 



IV. THE 3-POINT AND 4-POINT FUNCTIONS OF INFLATON PERTURBATIONS 



As an illustration of the formalism developed above, we will derive expressions for the tree-level 3-point and 4-point 
functions of the infiaton fluctuations. We shall work in the tensor gauge ([50)) . All of the relevant cubic vertices are 
presented in Appendix D in Eqs. (jl6H - ll64)) . Here we list the specific couplings that are relevant for computation of 
three and four point functions. 

From Eqs. (|161H164"|) we can easily extract the terms involving w^: 

.... - / d^-d. ^ { - ^) -'^ + ^^^.^) (1^^) + ^ ^ (^^) (^^) (93) 



Hw w + H -iP- H — H'^w^ H w[ 

Notice that the w'^w term can be integrated by parts to obtain terms that can be combined with the terms. The 
terms in Eqs. (I161H1641) of relevance for the tensor three-point function involving three gravitons are, 

S,,, = I d^xdtNa'l^ - \hin^^hr - iHhl^hl^kr - ^-^^^^hJ^hYh^ (94) 



There are also the cubic terms that involve two scalars and one graviton and those that involve one scalar and two 
gravitons, 

€.H ' n~"~p , T^'T' , , T"!^ 1 'T"T' ( didj 



e(7 + e — 2?7)e_ff^ ,ttiTT ^ uTT ('^^ uTt\ ^ f^'-^ij ' 



^r/'.5--H1i^^")-^(^)(^)r 



6 " »J y 6~ «j Va2 y 8 

Notice that the last term in ([M]) will contribute as 0{e'^) to the 4-point function 71]. For the calculation of the scalar 



The result I I93I1 can be compared to Eq. (3.8) of Ref. [13| , which, when expressed in our language, (f) — p H , ip ^ —{\/e/Z)w, 
V ~^ (~v^/3)('"' + (<= ~ v)Hw), V^x — ^ (f/3)io, reduces to. 



iMaldacena = J d^' xdt — ^{^— - j jww + —wid.w) l^ — W j + jW [^ — W j 1^ — 



J www I Maldacena 

J )J l \ Z D 

3 2 

V6 2 J V em 18H^ 6 / 6 

We extracted 1/9 in H93| l as opposed to 1/18 in l|93| l in order to compensate for the different definition of the Planck constant in 
Ref. [T3l . where 87rGjv = 1 (recall that here IGvrGjv = k-^ = 1). By comparing Eqs. I I93I I and II93I I we see that the leading order (0{e'^)) 
terms agree, while there are some disagreements in the subleading terms (recall that and V_,^0,^/\/e are both linear in slow roll 

parameters). Since it involves subdominant terms, this discrepancy is quantitatively unimportant. Nevertheless, it would be of interest 
to find the source of the disagreement. 
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four-point functions the terms involving w'^ and any of the auxihary fields {S, n, SN^ , p'-' , p^p} are also needed, 



wwlaux) 



(3 - e)Na'^H 



n 
N 



V 18 36/ 



e(e - ri) 



18 



18 



/ e 



,w~ ^ —Hww + I ]H w '^^w~ 

18 36/ 9 V2 9 18/ 



lis^ser" 



e (diW\(diW\ (OiOj .\(OiOj .\ ^/C^ 



18 



'did. 
9-"^^^lv 



-pw + 



{d^^Ni + fts, ^) . a,a, . X e{Nj + d,s) . 



/ftCj \ 



-w{diw) 



187Va2 V V2 / 9iVa2 

The purely quartic scalar terms can be read off from Eqs. (|170m72() in Appendix D 



wwww 



4 e^^^V^000 4 



1944 



972 



(97) 



(98) 



With the above interaction vertices one can, for example, compute the scalar contributions to the 3-point and 4-point 
functions, as we describe below. 

A. The 3-point function {uT') 

We can now use the diagrammatic rules described above to obtain the tree level 3-point function. The relevant 
diagrams are 



(99) 



where the vertex is determined by (|93p . Keeping the leading order (©(e^)) slow-roll terms in Eq. 

--^ ^ ft . \ /'diW\ /ft 



e 

36' 



and recalling Eqs. (15(11) . ([MHSSj) and ([S7|) we see that the first interaction term in (|100p contributes as 



+00 



iSl\'w{t,^)w{t,y)w{t,^)W^ = I / d^uiVdra^ — |(rj|rw;(t,x)u;(T,u)|l])_(r!|T?«(t,y)w;(T,u)|ri) 



(100) 



X (r2|Tw(t, z)w(t, u)|rj) + perms 
-(n|u;(T,u)w(t,x)|fi)_(r!|w;(r, u)u;(t, y)|rj) _(r!|u;(T, u)w;(i, z)|r!) -f perms | 

= jSndT a3_i_ I A-_ (r, u; i, x)9, A'ji^ (r, u; t, y)9, A'ji^ (r, u; i, z) + perms 



-A™ , (r, u; t, x)a^A'^: , (t, u; t, y)a^A'^: , (t, u; t, z) -I- perms 



(101) 
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which is just the retarded contribution. Including in the similar manner the other two leading order contributions 
from (|100p we finally get for the scalar three point function, 



{n\wit,^)w{t,y)w{t,z)\n) = yd3udT^|A!;i„(T,u;t,x 



1 d 
3 9m* 



A'^;+(r,u;t,x) 
d 1 



6a^ du^ 

drA%{T,u;t,y)drA'^^{T,u;t,z) 



perms 



(102) 



a.A-+(r,u;t,y)^^9.A-+(r,u;t,z) + ^^A-+(r,u;t,y)^A!:;+(r,u;t,z) 



d 

6a^ du'' 



' dw- 



perms 



4 



where perms denote the six permutations of {x, y,z}. When Eq. (|102l) is transformed to the spatial Fourier space, 
and inflationary mode functions are used for the scalar propagators, one can calculate the induced bispectrum, which 
in turn can be related to the cosmic background temperature fluctuations. As expected, the expression (jl02l) leads to 
a manifestly real result since A™_ = — (A!!i^)* , and coincides with that obtained in the operator formalism [l^ via 



{n\w{t,x)w{t,y)w{t,z)\n) = i / (fudT{n\[w{t,x.)w{t,y)w{t,z), 5i(T,u)]|rj) 



where Si is the interaction action in the tensor gauge and all fields are taken to be free Heisenberg fields. Thus, we 
have reproduced the results already obtained using the operator formalism and canonical quantization jl3l] . 

B. The 4-point Function {w'^) 

The quartic interaction terms in (j98p contribute diagrams of the form 



(103) 



+ 





However, the leading order contribution to the 4-point function comes from (j97p . corresponds to diagrams like 



(104) 



where the dashed line represents the "propagation" of auxiliary fields. As is evident from the form of the auxiliary 
propagators, such diagrams correspond to effective 4-point interactions. Since A-tzp = for all of the auxiliary fields, 
there can be no H — internal lines involving auxiliary fields. The leading order 0{e^) terms from Eq. (j97p . determining 
the vertices of the diagrams (|104p . read 



^wwlaux) 



d-^xdtNa- 



e n 
18N 

T 



tip -t- 



18 (3 - e)Na^H 



a 



e(iV, 
9 



Na 



diS) 

w[diW) 



(105) 
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Diagrams like (I104p contribute products of the form, 



ww{aux)'^ww{aux) 
2 



o 



81 



(wdiw) —r (wdiw) 



(106) 



324 



iwdrw)^^{wdjw) + 



324iJ 



d, 



n 



where O is an operator and we used the auxihary propagators (|74H76p . Note that the contributions from the two 
auxihary fields linking the terms h — V^S'/[(3 — e)a^H] cancel each other out. We see that the exchange of auxiliary 
fields induces an effective 4-point interaction given by the expression inside the curly brackets of (|106p . This effective 
interaction corresponds to the leading order terms of the quartic interaction terms derived in i72i . ,73j ■ 

The contribution to the four point function is obtained when O — > x)u)_|_(t, y)?i;_|_(t, z)w_|_(t, v). The result is, 



{n\w{t,x)w{t,y)w{t,z)w{t,v)\n) 



d^u dTNa^l—l^—A'f (r,u;t,x) — A!" (r, u; y)) — f ^— A!f (r,u;t 



z)^A^(r,u;t,v) 



648 

^2 



-Xa- 

NdT +- 



(r, u;t,x)^A!;;_(T,u;t,y) 



1 d d 



d 



d 



648i? 



A-„(T,u;t,x)_A-_(T,u;t,y) 



NdT 
d ( d 



Ndi 



( 



V2 du' du3 \NdT 

a 



r,u;t,z)^A™_(T,u;t,v) 



adu^ 



A-_(T,u;t,x 



d 



(+- 



adu^ 
+) + perms 



A-_(T,x;^,y) 



(107) 



where perms denote the 24 permutations of {x, y, z, v}. This leading order contribution is of the order 0{e^). There 
is also the tree level contribution coming from the graviton exchange. The leading order vertex in Eq. (|95p is of the 
order O(e^), and hence formally it also contributes at the order O(e^). Here we do not discuss this contribution any 
further, and refer the reader to Ref. 71] where the graviton exchange contribution to the 4-point function has been 
estimated. Finally, the terms in Eq. (|98p yield contributions that are of the order O(e^); since their contribution to 
the four point function is suppressed with respect to (|107p . we do not calculate them here. 



V. DISCUSSION 



In this paper we developed a path integral formulation for inflationary perturbations in single field inflation. Using a 
phase space formulation of the system as a starting point, we were able to bring the free action into a particularly simple 
form which allows for a straightforward calculation of the various propagators, recovering in the process standard 
results of linear gauge-invariant cosmological perturbation theory. The interaction terms can then be obtained without 
the necessity of first solving the energy and momentum constraints as has been done so far in the literature. A notable 
feature is the appearance of two types of auxiliary fields in the path integral beyond the physical propagating degrees 
of freedom: a set of commuting non-dynamical fields related to the existence of the constraints, as well as a set of 
anticommuting Faddeev-Popov ghost fields induced by the imposition of gauge conditions. The resulting path integral 
is independent of the choice of gauge, provided the asymptotic states are defined in terms of the linearized fields and 
the corresponding linearized gauge transformations. 

We then briefly described how to obtain N-point expectation values and commented on the meaning of two different 
gauges and their relation to cosmological observables in this formalism. They are the uniform field gauge and the 
tensor gauge which have been widely used in past considerations of non-Gaussianity but in the context of an interaction 
picture operator approach. Correlators can be expressed in a systematic expansion in diagrams in the in-in formalism. 
We found that, when quantum loops are taken into account, anticommuting ghost fields must be included in the 
computation. Furthermore, internal lines in diagrams should also include the commuting auxiliary fields and we 
demonstrated their role in the computation of 4-point functions, in which the leading order contributions indeed come 
from diagrams with internal lines involving auxiliary fields. This way of obtaining the effective 4-point interactions 
seems simpler than having to go through the solution of the constraints at second order and the substitution of these 
solutions back in the action [tI, [tII . 

So far, we have only considered standard single field infiation but the generalization to multi-field and more com- 
plicated models should be straightforward as long as the canonical formulation is known. We have also explicitly 
considered only tree diagrams in the tensor gauge, noting that the standard non-linear redefinitions (gauge trans- 
formations) can be used to obtain results for the curvature perturbation on comoving slices. Of course we could 
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have obtained this result by working directly in the uniform field gauge, at the price of more complicated interaction 
terms. Relating results in the two gauges would be particularly interesting if quantum loop corrections are taken into 
account and questions of renormalization arise. This is a regime where the role of the anticommuting ghosts would 
become crucial. Currently such questions are only addressed using the (essentially classical) AN formalism on long 
wavelengths. Other issues involve the appearance of infrared divergences, backreaction and a more rigorous definition 
of stochastic infiation. A path integral formulation might prove very useful for such considerations which will be the 
focus of future work. 
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Appendices 

In the following appendices we give details on the derivation of equations ([M)) and (|159p (Appendices A and B); 
the Poisson brackets needed to calculate (|5ip (Appendix C); and the general interaction hamiltonian as well an 
intermediate result needed for cubic and quartic vertices evaluated in section IV (Appendix D). 

Appendix A 

Let us consider the spatial Ricci scalar contribution the action ([TH3]). A spatial metric perturbation Squ — a^hi 
on a fiat FLRW background can be written as 



gij = a'^{6ij + hij) . 

Its matrix inverse g^^ = a^^g^-' is then 



(108) 



(109) 



where the indices are raised with the Kronecker delta, K^^ = 6^^hij, h^^ = S^^h/ . Eq. (|109|) can be also written in 
matrix notation as, 



g = I h + h h hhh 



(-l)"h" + C'(h"+^) 



(110) 



We will find the following notation useful: a numerical subscript n will indicate that terms of order n and higher in 
/ly are included. For example g>„ = (-l)"h" + 0(h"+i) etc. To calculate the determinant = [Det(I + h)]=^i/2, 

we can use 



[Det(I + h)]±V2 ^ exp±iTrln(I + h) = ^ iJ^T E 

71=0 ' \ m=l J 



(111) 



For example, in accordance with this notation, we have {h = Tr[h]) 



±iTr(h3) - i/iTr(h2) ± —h^ 
6 ^ ^ 8 ^ ^ 48 



(112) 



We will also be using 



Aiju = 6ikSji + 6u6jk — SijSki 



(113) 
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We can now write the spatial connection as 



(114) 



The corresponding Ricci scalar R = g^^Rij, with i?y = dk^^^ - djV'^,, + Tl^T^^ - r'j,^r^; is then, 

+^g'"'9''"{dkh„ii) {d,Kj + djh.n - dnhij) - ^g''"^''" (9fc/imj + djhkm - d^hkj) {^^Kl + dihin - a„/i,;)| . (115) 
Now making use of 

di~g'' = -ndihnmrg"'^ (116) 
and upon summing various identical terms, Eq. (jllSp can be written in the more symmetric form as, 



11 3 

'^{dlhjra){dkhn) - - ) (9; /lfo„) + {dlhij){dkhran) + -^{dihkl){djhmn) 



where (see Eq. (|112p ) 



(117) 



y/g = a^^/g 



h 



(118) 



Formula (|117p contains vertices to all orders generated by the ^/gR term, and it is written such that it is relatively 
easy to extract its contribution to an arbitrary n-point vertex. For example, the linear contribution in hij is 



i^/gRhin = a[didjhij ~ dfh] , 



(119) 



while the quadratic contribution reads. 



quad 



-{didjhij - V^ft.) - hki{2dkdihu - V^hki) + hijdidjh 

3 1 3 

-^{dkhki){dihii) - -{dih){dih) + {dih){djhji) + -{dihji){dihji) 

-^hdfh + ^hdidjhij - ^h^jdidikji + ^hjiW'^hji + (tot. der.) , 



(120) 



where in the last line we ignored several total spatial derivative terms (tot . der . ), which drop out when inserted into 
the (free part of the) action (P) . 



Appendix B 



We now present steps in the derivation of the free action P5|) and (I137P below. Let us first analyze the momentum 
part of the action ([T]-[3]) quadratic in perturbations (see Eq. p4l) ). which can be written in the form. 



where 



, . 2 

-h 09 H — 



tjAjjdki 



2H ' '■'2 N 



1 



1 ; 1, 3 

^ ^ 2H 2 N Na^H 



(121) 



(122) 
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It is now convenient to complete the square in Eq. (|121[) by the appropriate shifts of the momenta tt'-' and tt^ of the 
gravitational and scalar fields, respectively. From Eq. (|12ip we easily see that the shifted momenta are, 



With these we get for the momentum terms (|12ip 



(123) 
(124) 

(125) 



The linear shifts (I123P and (|124l) decouple the momenta from the inflaton and the metric perturbations. 

Next we observe from (|125p that the momenta shifts (|123l) and (|123p have introduced terms quadratic in the lapse 
perturbation n. Along with the original linear terms in the action (l23l) , they read 



d'^'xdt 



N 



n + n/„ 



where 



N 



(126) 



(127) 



Just like in the case of the momentum terms above, we now complete the square for the lapse. The contribution is, 



C(2) 



d-^xNdt 



N2 



N 



From the above manipulations we have the shift Ni also contributing to the quadratic action 



S 



(2) 



d xdt 



where 



and 



3V a6y; 



(5,0 ah 



2H 



P<4 2HV, 



(128) 



(129) 



(130) 



J = , (131) 



where in fact only the vector and scalar parts of Jij in (|13ip contribute. Analogously to the lapse, the shift contribution 
to the quadratic action can be written (up to boundary terms) in the form. 



S 



(2) 



~\ ■ 4H ~ / • \ 1 2-1 

h + 2VV^hj +2(f)H(p+ -^V^ih-V^h) + i^(t)V + 2HV^jip'^ |, (132) 



aN' 
3 . 1 



4(1 - a)V 



where 



1 



6(1 - a) 



J ~2a^NV^h 



(133) 
(134) 
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When all the terms from Eqs. (I125|) . (jl28l) . (I132[) and (p4)) are combined, we get the intermediate, complex expression 
for the action of quadratic perturbations in hij and (p (here we drop the contributions from the momenta and 
constraints) , 



c(2) 



(fxNdt- 



2H' 



5Tr[h • h] - -h^ 



hdfh — 2hdidihij + 2hijdidihji — hjidfhji 



Tr[h . h] + -h^ 



Tr[h • h] - /i^ 



8H 



Tr[h • h] - ^hh 



1 

'v 



2H' 

2 



8Tr[h • h] - 3h^ 



2Hh -~<j>ip+ -^{d^djh,, - dfh) - V^ifi 



(135) 



IT 12 



1 



2f4?- 4 df{h-dfh) 



(^<j)V + 2HV^^^ 



This action can be significantly simplified when expressed in terms of the transverse traceless tensor hfj'" defined in 
Eq. ([57H55| and the gauge invariant combination of scalar matter and scalar gravitational potentials. 



If = ip — z{h — V^/i) , 



(136) 



which is also known as the Sasaki-Mukhanov variable. An important step in the derivation of the quadratic action is 
a partial integration of the terms that contain one time derivative acting on perturbations ip and hij. Furthermore, 
since vectors are not sourced, their contribution to (jl35p cancels out, and the contributions from scalars and tensors 
separate. Following this procedure one eventually arrives at the quadratic action. 



Sf 



1 



iV2 



[^i^N, 



UP 



-J- 



3H- 



z z 



ij^? - (^) 



(137) 



where 



(138) 



and the shifts Lq, Li Lij and of the shifted momenta and constraints (pil - HH) in the first line of (|137p can be 
obtained from Eqs. (fT22HT23)) . (fHT)) and (fT33HT34)) and read 



N 
2a^V 
a^NV 



2a^Hh - V4,ip - a^V4,ip + a{^^^Jh^j - V^h) + 



aV 
6NH 



J -2Na^V^h 



12i72-'^ 24i/2v2 



AH 
Sii 
2 



1 • 2 

h + h ^Ln = — - — 

2H N Na^H 



OkLk - —n 



1 



N Na^H 



1-4- 



V 



-n+ -—(p = — + ^—z — V S , 

2 Vd, N N NaW 



2H 



V^S (139) 



(140) 



(141) 
(142) 



with 



VI ^^3^2iJP^ . , (V^ ^ 2HV 



(LI 
Va3 



V 



2H 



{dkdihki - V^h) 



(143) 
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We shall now argue that the action (|137|) is gauge invariant. To show this we need to consider how the metric and 
scalar transform under infinitesimal coordinate transformations, x'^ — >■ under which the metric transforms 

as, gf^u g^ii, — 2V(^^i,), = g^iy^". To evaluate the covariant derivative in the metric transformation need the 
Levi-Civita connection to the zcroth order in perturbations. The nonvanishing components are, 

J- 00 —^^1 ^ oi ~ ^ jO — Jvrroij, L — ^ ■ 
On the other hand, the metric components to linear order in perturbations are, 

g^,l^ = -2Nn; g^^^ = a' S,, = ; g^f=a^h,,, (144) 
where hij is decomposed as in Eqs. ([37H38)) . These then imply the following linear transformation rules, 

N 



^^^+^?o; h^h^2^^+^^o; hf^hf-^^f; h^T^h^T^ (145) 

TV a'^ N a'^ ■' 

where = S,f+di(, and di£,-f — 0. In particular, h — V^h h — V^h+{6H/N)(o- Now, from the transformations (|145p 
and equations (I136p . (|133H134|) . (|131l) and (|127l - |128|) it follows that hJJ^ , 'f', n, , S, and also the corresponding 
terms in the action (I137p . are gauge invariant. 

An important question is whether the whole action (jl37p is gauge invariant, also when the momentum terms are 
included. That the answer is yes can be argued as follows. When the Hamilton equations for the momenta, obtained 
by varying the action ([Ij, are solved in terms of the fields and inserted back into the action ([1]), one obtains the 
standard Einstein-Hilbert action. This then suggests that one can define the coordinate transformations for the 
momenta such that also the action ([1]) becomes covariant. When this program is carried through at linear order in 
coordinate transformations, one gets the transformation rules for the momenta which render the free action (jl37p 
gauge invariant. To see how this works in detail, note first that varying the action with respect to p^p and p^^ yields 
Pip = T^ip + Iip/2 = and p'^ = tt*^ + (/^ — 5ijI)/2 = 0, c/. Eqs. (|123H124|) . where I^, Itj and / = Sijltj are defined 
in Eqs. (|122p . These equations are just the Hamilton equations of the linearized theory and must be gauge invariant, 
implying that p^^ and p^p are gauge invariant. From the transformation rules 

+ ^V^e + ^Co ; ^ - + ^S,,V^C + ^^9.9,^0 + ^{H' + m,Co (146) 

and based on the Hamilton equations for the (unshifted) momenta tt^ = —1,^/2 and tt*^ = ^{^ij conclude 
that the momenta transform as. 

From this one can easily reconstruct how the scalar, vector and tensor components of tt*-' — {dij/3)TT + {didj — 
(5,;jV^/3)7r + d^iirj-^ + tt^'^, transform: 

o 

r^ij\TT . fij\TT. /^i\T . (i\T . ^ 

(tt-*) ^ {■k-') ; (tt ) ^ (tt ) + —ki 

^ ^ ^ + - ^T^TTTF^o ' '^^^ + 4^'^-7rT7TV% + T^[^^' + i?]eo. 

2Na^H Na''H NH '■ ■' 

The implications of these transformation rules for the cosmological perturbation theory are discussed in section IIIBI 
Next, the action (|137p is time reparametrization invariant. Thus choosing, for example, N = 1 ot N ^ a gives 
the action in physical (cosmological) or conformal time, respectively. One often defines the Mukhanov variable [5^ 
V = a(p. For example, when written in conformal time {N = a,t ^ rj) the lagrangian of the scalar part of the 
action (|137p will acquire an additional [a" /a]v'^ contribution (here a prime refers to a derivative w.r.t. conformal time 
r]): 

5(2) = / d'xdij- \v'^ - {d^v f + ^^vA . (148) 
J 2 [ az \ 
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Another convenient variable is 



z z 



in terms of which the scalar part of the quadratic action (|137p simplifies to, 



9(2) 



d,: 



6H 



gSign^, 



(149) 



(150) 



Appendix C 

We now give the expressions for the Poisson brackets of the constraints {Qq,(x), C^(y)}, which appear in the 
measure of the path integral (15^ and the constraints = C^'^ + C^^ are given in Eqs. ([55H55| and (|i5] - |T7l) . with 
the free and interaction hamiltonian given in Eqs. (1341) and (I159p . respectively. Before we proceed, we first present 
the Poisson brackets between the fields and the constraints, 



{/i„(x,i),Co(y,t),} 



'2H{g ^/^)|(5.y + [-h^j + h6,j - + 7r%] + \hh^ - huhij - Air'^hij + nh,, + n'^^hkAj] 
+ [-2h,y^hkj + hkyh,,]'j6{x - y) (151) 



}^(x 



{/i„(x,<),G(y,i)} ^ 

Mx,t),Co(y,t),} - -^(r'/')(l + 7r^)<5(x-y) 
{^(x,t),a(y,t),} : 



a 



(152) 
(153) 
(154) 



We are now ready to consider different gauges, 
a) "Tensor Gauge": Qo ^ h , Q,^ dj (h,,j - 



{Qo(x),Co(y)} 
{Qo(x),Q(y)} 
{Q,(x),Co(y)} 

{Q,(x),C,(y)} 



2H{g'^) 3 + (2/i + tt) + {-hijhij + - n^'hij + Trh) + {~2hikn''^hu + n^^h.^h) (5(x - y) 



a 



2 r 



S, 



(155) 



hiihij 



-^{%V2 + iar^J - ~9'\2d(,hi)k - d,hu)iy,t) or + lg'''{2dihik - duh){y,t) af }5(x - y) . 



Since the tensor gauge conditions Qa — set the spatial scalars and vectors to zero, h = = h = hf , from the 
relations (|155p one easily obtains the ghost operators flap — {Qa, C^}|q„=o of Eq. ([STj) : 



~2H{g-"^f^ [i + P + L- hl^hl^ - {p^^ + U,)hl^ - 2hJ^{p^' + Lki)hl^\ <5(x - y) 
-\df5{^-y) 

-2H{r-r[-hJ^-2(p^^+U, ^f{p + L))Hp + L)hl^-{hf,^hl^-^fhirhir 



(156) 



{p'^^ + L,,)hlfhJ^}{y,t)d^5{^-y) 

^{%v^ + iara; - {g^r^{d.hlf + d^h^J - d,hfniy,t) ar}5(x - y) , 
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where L = LijSij is defined in (|14ip . where h = = h = hf and ip —[P,f,/{6a^H)]w are to be exacted, resulting in 
the shift functions given in Eqs. (|164H168() below. Finally, from Eqs. (|109p and (|112p we infer 



- hf;' + hThfi' 



,TTiTT-lTT 
'''il '''Ik l^kj 



^ hTT i^TT iTT iTT 

-n^j nji rill, ''ki 
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i,TT i,TT iTT i,TT 
"ij "ij "fci "fei 



7f) 



b) "Uniform Field Gauge" : Qo = V , = dj [htj - ■ 

In this gauge the relevant Poisson brackets follow from Eqs. (jl53H154")) . 



{Qo(x,t),Co(y,t)} 



1 



g 2(l + 7r^)5(x-y) 



{Qo(x,t),a(y,t)} = g'J(a^.^)^(x-y) 



(157) 



and the other Poisson brackets {Qi(x, t), C;3(y, i)} are identical to the corresponding tensor gauge expressions in 
Eq. ()155|) . In this gauge Lp = = h = hf , such that hij = {Sij/3)h + hf^', and the corresponding ghosts operators 
are, 



^00 



_2iJ(,g-^)"fe{-/,^.^-2(p- - ^fip + L)) + {lh + p + L)hJ/^-{hrhf;'-^fhffhl^) 



(158) 



4/i 

T 



rl 



(p + i)/i+(/'+LfcO^' 



TT 



}(y,05J5( 



X- yj 



-{g^^T\d.,hf^ + - dkhff'^Y, t) - i [(.g^"'=)"f«afe/.,](y, t) 9f }j(x - y) , 

where the superscript in (5^^^^)"^^ and (^*-')"fs signify that only h and /i^^ contribute to /ly , i.e. /i^j = {Sij /3)h + h. 



TT 
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Appendix D: Interaction Hamiltonian, cubic and quartic action 



The interaction hamiltonian containing cubic and higher order interactions as inferred from Eqs. ([2]) and (jll7|) is 
of the general form: 



3 1 1 1 1 

^{rh>3 - {rh>2h + h,j{g-^)>i-A,jkihki 



kl 



I 1 °° T/(") 



„l ^ 



n^2 



n— 3 



"P / 1 

■aN ((g^)>i ((5'^)>i<5'=' + (g"')>i) + 5^ (5*^')>i(.9")>i) {dAhj - d^d^ 

11 3 

--{dihj,n){dkhin) - -{dnhij){dihkm) + {dihij){dkh„in) + -^{dihki){djh„ 



aN 



(159) 



where gg"^ and (g±i/2)(>n) 

are given in Eqs. (|109p and (|112p . respectively. The cubic part of the interaction 

action (IMl) is of the form, 



(160) 



where Cq^^ and c|^'' are the quadratic part of Eqs. (|46ll47p and (p4| . When written in tensor gauge (|80)) . the cubic 
part of the interaction hamiltonian density (|159p reads 



Hcubic = ^a3<j2ij2 



TT 



,TT 



1 

4' 



18a2 y 

while the cubic constraint contributions read, 



(3-7?)6H^ ^^^r e^/^l^^ 
6 162 



N.c'f^ = {N,+L,)a- 



TT 



(161) 



(162) 



2H\{p^ + L 



(2eH 



-w 



I 3a2 



3a2 



(163) 
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The shift functions (IT39HI431) are, 



^0 

Lip 
Lij 
L 



Ne 



6 {3-e)a^H 
a^Ne a . 



6 V2 

w rjw n 
"m ^ ~6~ ~ iV ^ (3 - e)Na^H 



'J 12H 



3n 



6 



2Na?H N 2(3 - e)Na'^H 



e . e an V^S" 
w H — w =- + — - — 

6H 2 N {3-e)Na^H 



(164) 
(165) 
(166) 

(167) 
(168) 



Finally, for the calculation of four point functions, the quartic contribution from the interaction action (|24l) is needed. 



5, 



quartic 



d'^xdt 



{ 



^quartic H~ f^Ci 



(3) 




(169) 



where Cq^"* and Cl''' are the cubic parts of Eqs. (|46ll47p . In tensor gauge they are. 



.(3) 



nC, 



(3) 



l^cubic (170) 

while the quartic part of the hamiltonian (|159p becomes. 



K 



quartic 



uTTuTTi^TTuTT . uTT uTT uTT i,TT 
'^ij ''■jl ''■Ik 'T'ki + '^''ij "-y "fci "fei 



SijSknihip hp„ + 5ijhf.p hp^Sin + hip h^^ SkmSin + Sijh^^h^i + h^j SkmKi 



2eH' 
e 

18 



~'''kl 'T'kl "ij 



V a / V a 



1944 



72 



-whi/ hi' hi. 



{3-e)H\\h^;^hl^hlfhr - ^^^^ 



(172) 



These formulae are used in section IV to construct some of the cubic and quartic vertices of the theory. 
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